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ABSTRACT 



Context. Observations of molecular spectral lines provide information on the gas kinematics and chemistry of star-forming regions. 
Aims. We attempt to achieve a better understanding of the gas distribution and velocity field around the deeply embedded Class 
protostar SMM 3 in the Orion B9 star-forming region. 

Methods. Using the APEX 12-m telescope, we mapped the line emission from the 7 = 2-1 rotational transition of two CO isotopo- 
logues, 13 CO and C 18 0, over a4'x4' region around Orion B9/SMM 3. 

Results. Both the 13 CO and C 18 lines exhibit two well separated velocity components at about 1.3 and 8.7 km s -1 . The emission 
of both CO isotopologues is more widely distributed than the submillimetre dust continuum emission as probed by LAB OCA. The 
LAB OCA 870-yt/m peak position of SMM 3 is devoid of strong CO isotopologue emission, which is consistent with our earlier de- 
tection of strong CO depletion in the source. No signatures of a large-scale outflow were found towards SMM 3. The 13 CO and C 18 
emission seen at ~ 1.3 km s _1 is concentrated into a single clump-like feature at the eastern part of the map. The peak H 2 column 
density towards a C 18 maximum of the low- velocity component is estimated to be ~ 10 22 cm -2 . A velocity gradient was found across 
both the 13 CO and C 18 maps. Interestingly, SMM 3 lies on the border of this velocity gradient. 

Conclusions. The 13 CO and C 18 emission at ~ 1.3 km s _1 is likely to originate from the "low- velocity part" of Orion B. Our analysis 
suggests that it contains high density gas (~ 10 22 H 2 molecules per cm 2 ), which conforms to our earlier detection of deuterated species 
at similarly low radial velocities. Higher-resolution observations would be needed to clarify the outflow activity of SMM 3. The sharp 
velocity gradient in the region might represent a shock front resulting from the feedback from the nearby expanding H n region NGC 
2024. The formation of SMM 3, and possibly of the other members of Orion B9, might have been triggered by this feedback. 
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1. Introduction 

The protostellar phase of low-mass star formation be- 
gins when a starless (prestellar) core collapses, and, af- 
ter a hypothesised short-lived first-hydrostatic core stage 
(ILarson 1969I [Masunag a et al. 1998| ), a stellar embryo forms 
in its centre (the so-calle d second hydrostatic core; e.g., 
Masunaga & Inutsuka 2000 ). The dense cores harbouring the 
youngest protostars are known as the Class objects 
(I Andre et al. 1993 1 2000). In these objects, most of the sys- 
tem's mass resides in the dense envelope, i.e., M env » 
M*, where M+ is the mass of the central protostar. For 
this reason, Class objects, or at least the youngest of 
them, are expected to still represent the initial physical con- 
ditions prevailing at the time of collapse phase. Class ob- 
jects are characterised by accretion-powered jets and mole- 
cular outflows, which can be very powerful and highly col- 
limated (e.g., |Bontemps et al. 1996"! IGueth & Guilloteau 19991 
Arce & Sargent 2005| 2006; ILee et al. 20071) . The statistical life- 
time of the Class stage is estimated to be ~ 1 x 10 5 yr 



Send offprint requests to: O. Miettinen 

* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) under programme 088.F-9311A. 
APEX is a collaboration between the Max-Planck-Institut fur 
Radioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory. 



(Eva ns et al. 20091 Enoch et al. 2009b , but the exact duration of 
this embedded phase of evolution can be highly dependent on 
the initial/environmental conditions (e.g., |Vorobyov 20 it)} . 

The target source of the present study is the Class protostel- 
lar core SMM 3 in the Orion B9 star-forming region, which was 
originally discovered by Miettinen et al. (2009; Paper I) through 
LABOCA 870-yum dust continuum mapping of the region. SMM 
3 is a strong submm emitting dust core (^870 - 2.5 Jy) that is as- 
sociated with a weak Spitzer 24-pm point source (S24 - 5 mJy), 
and a 3.6 Jy point source at 70 pm. Using the EfFelsberg 100- 
m telescope NH3 observations, Miettinen et al. (2010; Paper II) 
derived the gas kinetic temperature of Ty± n = 11.3 ±0.8 K in 
SMM 3. Using this temperature, the core mass was determined 
to be 7.8 ± 1.6 M©, and its volume-averaged H2 number den- 
sity was estimated to be 1.1 ± 0.2 x 10 5 cm" 3 . In the SABOCA 
350-yum mapping of Orion B9 by Miettinen et al. (2012; Paper 
III), SMM 3 was found to be by far the strongest source in the 
mapped area (S350 - 5.4 Jy). We also found that it contains two 
subfragments, or condensations (we called SMM 3b and 3c), ly- 
ing about 36" - 51" in projection from the central protostar. 
These correspond to 0.08-0.11 pc or ~ 1.7 - 2.3 x 10 4 AU at 
d - 450 pcQ. Because the thermal Jeans length of the core is 
A] = 0.07 pc, we suggested that the core fragmentation into con- 



1 In this paper, we adopt a distance o f 450 pc to the Orion giant 
molecular cloud (Genzel & Stutzki 1989 ). The actual distance may be 



1 



Miettinen: 13 C0 and C 18 mapping of Ori B9-SMM 3 



densations can be explained by thermal Jeans instability. Using 
the 350/870-yum flux density ratio, we determined the dust tem- 
perature of the core to be T^ ust = 10.8^ K, which is very close 
to 7km within the error bars. The revised spectral energy distri- 
bution (SED) of the core yielded a very low dust temperature of 
8 K, and a bolometric luminosity of Lb i = 1.2 + 0.1 L©. The lat- 
ter is very close to the median luminosity of protostars in nearby 
star-forming regions, i.e., L me( j = 1.5^4 L© (Enoch et al. 2009; 
Offn er & McKee 201TT) . In Paper III, we also studied the che- 
mistry of SMM 3. We derived a large CO depletion factor of 
/d(CO) = 10.8 + 2.2, and a high level of deuterium fractiona- 
tion, i.e., a N 2 D + /N 2 H + column density ratio of 0.338 + 0.092. 
In Fig. CD we show the LABOCA 870-//m, SABOCA 350-yum, 
and Spitzer 4.5/24-yum images towards SMM 3. In Table [U we 
provide an overview of the physical and chemical properties of 
SMM 3 derived in our previous papers. 

In this paper, we discuss the results of our 13 CO and C 18 
mapping observations of the environment of SMM 3. We ana- 
lyse the structure of the mapped region as traced by emission 
from the / = 2 - 1 rotational transition of the above CO iso- 
topologues. The rest of the present paper is organised as fol- 
lows. Observations and data reduction are described in Sect. 2. 
Mapping results and analysis are presented in Sect. 3. In Sect. 4, 
we discuss our results, and in Sect. 5, we summarise and con- 
clude the paper. 
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2. Observations and data reduction 

The observations presented in this paper were made on 13 
November 2011 using the APEX 12-m telescope located at 
Llano de Chajnantor in the Atacama desert of Chile. The tele- 
scope and its performance are described in the paper by Glisten 
et al. (2006). An area of 4' x 4' (0.52 pc x 0.52 pc at d = 450 pc) 
was simultaneously mapped in the / = 2 - 1 rotational lines 
of 13 CO and C 18 using the total power on-the-fly mode to- 
wards SMM 3 centred on the coordinates a 2 oooo = 05 h 42 m 45 s .8, 
£20000 = -01° 16' 1370 [(+775, +370) offset from the SABOCA 
peak position of SMM3]. At the 13 CO(2 - 1) and C 18 0(2 - 1) 
line frequencies, 220 398.70056 and 219 560.357 MH2B, respec- 
tively, the telescope beam size is about 2873 (HPBW). The target 
area was scanned alternately in right ascension and declination, 
i.e., in zigzags to ensure minimal striping artefacts in the final 
data cubes. Both the stepsize between the subscans and the an- 
gular separation between two successive dumps was 974, i.e., 
about 1/3 times the beam HPBW ensuring Nyquist sampling. 
We note that the readout spacing 1/3 x HPBW should not be ex- 
ceeded to avoid beam smearing. The integration time per dump 
and per pixel was 1 s. 

As a frontend, we used the APEX-1 receiver of the Swedish 
Heterodyne Facility Instrument (SHeFI; Belitsky et al. 2007; 
IVassilev et al. 2008al b). The backend was the RPG extended 
bandwidth Fast F ourier Transfrom Spectrometer (XFFTS; cf. 
Klei n et al. 2012b with an instantaneous bandwidth of 2.5 GHz 
and 32768 spectral channels. The resulting channel separation, 
76.3 kHz, corresponds to about 0.1 km s" 1 at 220 GHz. 

The telescope pointing accuracy was checked by CO(2 - 1) 
cross maps of the variable star RAFGL865 (VI 259 Ori), and 



Spitzer I R AC MIPS 
4.5 & 24 |jm 



somewhat smaller as, for example, Menten et al. (2007) determined a 
trigonometric parallax distance of 414 + 7 pc to the Orion Nebula. 

2 The 13 CO(2 - 1) frequency was taken from Cazzoli et al. (2004), 
and it refers to the strongest hyperfine component F = 5/2-3/2. The 
C 18 0(2- 1) frequency was adopted fro m the JPL spectroscopic database 
at http : //spec . jpl . nasa . gov/ (Picke tt et al. 19981 ). 




Right ascension (J2000) 



Fig. 1. LABOCA 870-yum (top), SABOCA 350-yum (middle), 
and a Spitzer IRAC/MIPS two-colour composite image (bottom', 
4.5 fim in green and 24 fim in red) of the Class protostellar 
core SMM 3 in Orion B9. The LABOCA and Spitzer images are 
shown with linear scaling, while the SABOCA image is shown 
with a square-root scaling to improve the contrast between bright 
and faint features. The LABOCA contours, plotted in white, go 
from 0.1 (~ 3.3<t) to 1.0 Jy beam -1 in steps of 0.1 Jy beam -1 . 
The SABOCA contour levels, plotted in green, start at 3cr and 
are 0.18 Jy beam" 1 x [1, 2, 4, 6, 8, 10, 12, 14, 16]. In the bot- 
tom panel, the first SABOCA contour, i.e., the 3cr emission level, 
is plotted in white to guide the eye, and the white cross indicates 
the SABOCA peak position of SMM 3. The small subcondensa- 
tions SMM 3b and 3c discovered in Paper III are labeled in the 
middle panel. The green plus sign shows the target position of 
our previous molecular-line observations (i.e., the submm peak 
position of the LABOCA map before adjusting the pointing; see 
Paper III for details). A scale bar indicating the 0.05 pc pro- 
jected length is shown in the bottom left of the top panel, with 
the assumption of a 450 pc line-of-sight distance. The effective 
LABOCA and SABOCA beams, ~ 20" and 1076, are shown in 
the lower right corners of the corresponding panels. 



was found to be consistent within < 4". The focus was checked 
by measurements on Jupiter. Calibration was made by means of 
the chopper- wheel technique and the output intensity scale given 
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Table 1. Summary of the properties of SMM 3. 



Parameter 



Value 



SMM 3 

^2000.0° 
#2000.0^ 
VLSR^ 

rd 
kin 

re 
dust 
y^SED,cold 

dUSt d 
<T NT 

o~m/c s d 
M f 

N(R 2 y 

<n(H 2 )>/ 

J bol — -^cold -^warr 
^submm/-^bol 

/d(CO) 
[N 2 D + ]/[N 2 H + ] 
SMM 3b 

^2000.0 a 
^2000.0° 

N(u 2 y 

SMM 3c 

^2000.0 a 
^2000.0° 

N(n 2 y 



05 h 42 m 45 s .3 
-01°16'16" 
8.68 + 0.06 km s" 1 
13'/ 3 (0.03 pc) 
11.3 + 0.8 K 
10.8^ K 
8.0 K 
0.14 + 0.003 km s" 1 
0.7 + 0.04 
7.8 ±1.6 Mo/2. 1 ±0.8 Mo 
0.5 + 0.1 

8.4 + 1.1 x 10 22 /1.0 + 0.3 x 10 23 cm" 2 
1.1 + 0.2 x 10 5 /4.0 + 1.5 x 10 5 cm" 3 
(0.3 + 0.1) + (0.9 + 0.1) = 1.2 + 0.1 L Q 
0.1 
10.8 + 2.2 
0.338 + 0.092 

05 h 42 m 47.6 
-01°16'24" 
0.7 + 0.2 x 10 22 cm" 2 

05 h 42 m 48 s .6 
-01°16'32" 
0.7 + 0.2 x 10 22 cm" 2 



Notes. {a) SABOCA 350-//m peak position. w The LSR velocity derived 
from optically thin C 17 0(2 - 1) line. (c) Effective radius of the "main" 
core as determined from the SABOCA 350-//m map. (J) Derived from 
NH 3 data. <x NT and c s are, respectively, the one dimensional non-thermal 
velocity dispersion and the isothermal sound speed. (e) Computed from 
the 350-to-870 fim flux density ratio. (jf) The first value refers to the 
LAB OCA 870-yum core, and the second one to the "main" core detected 
at 350 fim. ig) Virial parameter defined by Qf vir = M vir /M. (/l) L submm is 
the submm luminosity derived by integrating the SED longward of 350 
/umS l) Calculated by making the assumption that T dust equals the 7\ in 
derived for the "main" core. 



by the system is T*, which represents the antenna temperature 
corrected for the atmospheric attenuation. The amount of pre- 
cipitable water vapour (PWV) was in the range 1.28 - 1.48 mm, 
and the single- sideband system temperature was around 150 K 
(in units). The observed intensities were converted to the 
main-beam brightness temperature scale by Tmb = T*/tj M b, 
where tjmb = 0.75 is the main-beam efficiency at the frequen- 
cies used. The absolute calibration uncertainty is estimated to be 
about 10%. 

The spectra were reduced and the maps were produced us- 
ing the CLASS90 and GREG programmes of the GILDAS soft- 
ware packag^E The individual spectra were Hanning- smoothed 
to improve the signal-to-noise ratio of the data. A third-order 
polynomial was applied to correct the baseline in the spectra. 
The resulting Icr rms noise level of the average spectra are about 
30 mK (in Tmb) at the smoothed resolution (16384 channels). 
The data were convolved with a Gaussian of 1/3 times the beam 
HPB W, and therefore the effective angular resolution of the final 
data cubes is about 30". The average Icr rms noise level of the 
completed maps ranges from 0.20 to 0.28 K per 0.2 km s" 1 chan- 



nel on a Tmb scale. The constructed data cubes were exported in 
FITS format for further processing in IDL. 



3. Mapping results and analysis 

3.1. The average spectra 

The average 13 CO(2 - 1) and C 18 0(2 - 1) spectra are shown 
in Fig. [2] Both lines exhibit two well- separated velocity com- 
ponents: one near the systemic velocity of about 8.7 km s" 1 , 
and the other at 1.3-1.4 km s" 1 . It is not surprising that we 
see these lower- velocity components in the lines of CO isotopo- 
logues. The additional lines at comparable radial velocities of 
1.3-1.9 km s" 1 were already detected in the lines of N2H + (1 - 0) 
and N 2 D + (2 - 1) in Paper I, NH 3 (1, 1) and (2, 2) in Paper II, 
and C 17 0(2 - 1), H 13 CO + (4 - 3), DCO + (4 - 3), N 2 H + (3 - 2), 
and N 2 D + (3 - 2) in Paper III towards other cores in Orion B9. 
Therefore, detection of lower-velocity line emission from CO 
isotopologues was expected. We note that the average 13 CO line 
near the systemic velocity of SMM 3 appears to show a blue 
asymmetric profile with blue peak being stronger than the red 
peak. The central dip also appears to be near the radial veloc- 
ity derived from optically thin C 17 0(2 - 1) line in Paper III. 
Despite of these characteristics, the double-peaked li ne profile is 
not caused by infall motions (e.g., Myers et al. 1996| ); it results 
from averaging over the entire mapped area, where two separate 
velocity components at about 8.5 and 9.5 km s" 1 are seen. This 
issue will be further discussed in Sect. 3.3. The average C 18 
line profile at the systemic velocity is nearly Gaussian, which 
suggests that the line is likely to be optically thin. However, a 
hint of the two nearby velocity components is also visible in 
the average C 18 line; the line exhibits a small "knee" at ~ 9.5 
km s" 1 . 

It can also be seen from Fig. [2] that some of the observation 
OFF positions had 13 CO emission in the velocity regime bet- 
ween about 2.4 and 4.8 km s" 1 , and between about 10.3 and 1 1.3 
km s" 1 , and C 18 emission in the range 2.0-3.3 km s" 1 . These 
velocity regimes show up as artificial absorption features in the 
average spectra. Given the ubiquitous nature of multiple velocity 
components along the line of sight towards Orion B9, finding an 
emission free OFF position from this region can be difficult. 

The main purpose of examining the average spectra is to 
determine the velocity range of the detected emission. This is 
needed to construct the line-emission maps, as will be described 
in the next section. 



- Average spectra ft 

















3 Grenoble Image and Line Data Analysis Software is pro- 
vided and actively developed by IRAM, and is available at 



Fig. 2. Hanning- smoothed spatially averaged 13 CO(2 - 1) and 
C 18 0(2 - 1) spectra across the mapped field. The C 18 0(2 - 1) 
spectrum is offset by -1.45 K from zero baseline for reasons of 
clarity. The vertical red dashed line indicates the systemic veloc- 
ity of SMM 3 derived from C 17 0(2 - 1) in Paper III. 



|http : //www . iram . f r/IRAMFR/GILDAS 
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3.2. Moment maps 

To display the intensity and kinematic structure of the 13 CO 
and C 18 line emission, we constructed the moment maps by 
integrating the lines over the following LSR velocity ranges: 
[6.09, 10.36] km s" 1 and [7.77, 10.36] km s" 1 for the 13 CO and 
C 18 lines of the main velocity component, and [-0.07, 2.52] 
km s" 1 and [0.14, 2.10] km s" 1 for the 13 CO and C 18 lines 
of the lower- velocity component (at ~ 1.3 km s" 1 ). These line 
windows were selected so that the artificial absorption features 
discussed above are avoided. The threshold used for the moment 
maps was chosen to be 2 times the rms noise, i.e., 2<x. 

The zeroth, first, and second moment maps (i.e., the images 
of integrated intensity, intensity-weighted central velocity, and 
intensity-weighted FWHM linewidth) of the 13 CO and C 18 
emission of the main velocity component are shown in Fig. [3] 
while those of the lower- velocity component are shown in Fig.|4l 

The map of 13 CO integrated intensity (upper left panel of 
Fig. 0) shows that there is a local emission minimum close 
to the SMM 3 protostellar position. This conforms to the high 
level of CO depletion derived from C 17 data. The 13 CO emis- 
sion around SMM 3 is rather extended, which is not surpris- 
ing because it arises from the lower-density gas. In contrast, the 
LAB OCA dust continuum emission shows only the densest part 
of the region, i.e., SMM 3. The 13 CO emission appears to be 
strongest at about 2' south of SMM 3, and from there it extends 
to the northwest part of the map. From the lower left panel of 
Fig- El it can be seen that the C 18 emission follows quite well 
the morphology of the 13 CO emission. There is a hint of an elon- 
gated filament-type feature along the NE to the SW direction 
with relatively strong emission. There appears to be a few C 18 
maxima at the eastern part of SMM 3, one lying at the eastern 
tip of the 3.3cr LABOCA contour, and the other about 24" from 
the central protostar. 

In the case of the lower- velocity component, both the 13 CO 
and C 18 emission are less extended, and are instead concen- 
trated into a single clump-like feature at the eastern part of the 
mapped region (left panels of Fig.©. There is a 13 CO extension 
to the south of SMM 3 in projection, which is not seen in C 18 0. 

Interestingly, as can be seen from the first-order moment map 
of 13 CO (top middle panel of Fig. [3]), there is a fairly sharp bor- 
der between the two velocity fields, and SMM 3 appears to lie 
exactly between them, i.e., at the border of the velocity gradient. 
There is also a hint of increasing 13 CO linewidth across this bor- 
der as shown in the top right panel of Fig.(3j The radial- velocity 
structure of C 18 emission is quite similar to that of 13 CO, but 
the C 18 linewidths show less obvious spatial trend (lower mid- 
dle and right panels of Fig. [3]). 

Another interesting feature concerning the radial-velocity 
distribution is that also the lower- velocity component of C 18 
shows a somewhat similar gradient across the map (Fig.|4l lower 
middle panel). Implications of the velocity gradient will be dis- 
cussed further in Sect. 4.3. 

3.3. Channel maps and PV plots 

Velocity channel maps of the 13 CO and C 18 emission for the 
main velocity component are plotted in Figs. [5] and respec- 
tively. As shown in these maps, extended 13 CO emission can be 
seen across the velocity range 7.86 < vlsr < 10.15 km s" 1 , 
while that of C 18 extends over a somewhat narrower range of 
velocities, 8.32 < v L sr < 9.78 km s" 1 . 

Figure [7] shows a position-velocity (PV) diagram for 
13 CO(2 - 1) emission with the direction of the cut shown in 



Fig. [3] The slice is taken through the entire mapped area, ex- 
tending from the northwest corner to southeast corner, and its 
position angle, measured east of north, is P.A.= 135°. As shown 
in the top middle panel of Fig. O the PV slice goes across the 
border of the velocity gradient. The presence of two velocity 
components at ~ 8.5 and ~ 9.5 km s" 1 are clearly visible in the 
PV plot. These velocities correspond to the NW and SE parts of 
the mapped area, respectively [the PV diagram of C 18 emission 
(not shown) is essentially similar] . 




5 h 42 m 50* 40* 
R.A. 



Fig. 5. Velocity channel map across 6-10.56 km s" 1 of the 
13 CO(2- 1) line towards SMM 3. The velocity of each channel is 
shown in the top-left corner of each panel. The solid contours go 
from 0.67 to 4.67 K in steps of 0.67 K (Tmb), while the dashed 
contours mirror negative values due to the baseline-fitting prob- 
lems. The central cross on each panel indicates the map centre. 



3.4. Spectra, line parameters, and column densities at the 
selected positions 

We extracted individual 13 CO and C 18 spectra from the map 
centre, i.e., near the LABOCA 870-yum peak of SMM 3. 
Moreover, to investigate the properties of the lower-velocity 
component, we extracted the 13 CO and C 18 spectra from a 
peak position of C 18 integrated intensity lying east of SMM 
3 in projection (see the left panels of Fig.|4|. The extracted spec- 
tra are shown in Fig. [8] The C 18 line towards SMM 3 shows 
two nearby velocity components. There is a hint of that also 
in the 13 CO spectrum. Moreover, the 13 CO line appears to be 
blueshifted with respect to the C 18 line, suggesting that the 
lines originate in two different gas layers along the line of sight. 
In the case of the lower velocity-component spectra (lower panel 
of Fig. [8]), the 13 CO line shows two velocity components around 
the systemic velocity of 1.38 km s" 1 derived from C 18 0(2 - 1). 

The line parameters of the extracted spectra are listed in 
Table [2] In this table, we give the offset of the target position 
from the map centre (in arcsec), radial LSR velocity (vlsr), 
FWHM linewidth (Av), peak intensity (Tmb), integrated inten- 
sity ( J 7mb<1v), peak line optical-thickness (tq), and excitation 
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Fig. 3. The colour scale in the panels from left to right shows, respectively, the Oth (integrated intensity in Tmb), 1st (intensity- 
weighted centroid velocity structure), and 2nd moment (intensity-weighted linewidth structure) map of the 8.7-km s" 1 component 
13 CO(2 - 1) (top row) and C 18 0(2 - 1) emission (bottom row). The black contours in the upper panels, from left to right, go as 
follows: from 0.7 (3.2cr) to 7.0 K km s" 1 (32cr) in steps of 0.7 K km s" 1 ; from 7 to 9.5 km s" 1 in steps of 0.5 km s" 1 ; and from 1 
to 4 km s" 1 in steps of 1 km s" 1 . The corresponding contours in the lower panels are from 0.4 (2<x) to 3.2 K km s" 1 (16<x) in steps 
of 0.4 K km s" 1 ; from 7 to 9.5 km s" 1 in steps of 0.5 km s" 1 ; and from 0.5 to 2 km s" 1 in steps of 0.5 km s" 1 . Superimposed on the 
maps are the LAB OCA 870-yum contours in white as in Fig. [TJ The yellow solid line on the top left and middle panels shows the 
location of the PV slice extracted for Fig.[7J The effective beam size of 30" is shown in the upper left corner of the upper left panel. 



temperature (r ex ). The values of vlsr and Av of the 13 CO(2 - 1) 
lines were derived through fitting the hyperfine structure of the 
line (see lCazzoli et al. 2004b . while those of C 18 0(2 - 1) were 
obtained by fitting the lines with a single Gaussian. The two 
nearby velocity components were fitted simultaneously but for 
SMM 3 we only give the parameters of the components overlaid 
with green lines in Fig. [8] In contrast, in the case of the C 18 O 
peak of the lower velocity-component, we list the 13 CO parame- 
ters of the component slightly redshifted from the C 18 line ve- 
locity. The integrated line intensities listed in Col. (7) of Tabled 
were computed over the velocity range given in square brack- 
ets in the corresponding column. This way, we were able to take 
the non-Gaussian features of some of the lines into account. The 
quoted uncertainties in vlsr and Av are formal Icr fitting errors, 
while those in Tmb and J Tmb^v were estimated by summing in 
quadrature the fitting error and the 10% calibration uncertainty. 

The values of To and T QX for the lines towards the C 18 peak 
of the lower velocity-component were derived as follows. By 
making the assumption that the 13 CO and C 18 emission arise 
from the same gaqj and the two transitions have the same beam 



4 The LSR velocities of the two transitions are slightly different from 
each other [see Col. (4) of Table |2), so they may not be tracing exactly 
the same gas layers. 



filling factor and excitation temperature^, we can numerically 
estimate the line optical thicknesses from the ratio of the line 
peak intensities as (e.g., M yers et al. 1983||Ladd et al~1 998 ) 



r MB ( 13 CO) ^ 1 - 

r MB (c 18 o) * 1-e-is 



(i) 



The peak intensities are used here rather than the integrated in- 
tensities because multiple velocity components could "contam- 
inate" the latter values. To calculate the transition optical thick- 
nesses, we used the CO isotopologue abundance ratio of 



[ 13 CO] = [ 13 C] [ 16 Q] 1 
[C 18 0] [ 12 C][ 18 0] 60 



8.3. 



(2) 



The adopted carbon- and oxygen-isotopic ratios are the same as 
those used in Paper III for a proper comparison to that work 
(see references thereinfl Consequently, the optical thickness ra- 
tio T13 /ri8 was assumed to be equal to 8.3. 



5 We note that the observed 13 CO and C 18 transitions have similar 
frequencies. Therefore, the frequency-dependent main-beam efficiency 
(?7mbX and the telescope beam HPBW are also almost identical for the 
two transitions. 

6 We note that larger rat ios of [ 12 C]/[ 13 C] = 77 and [ 16 0]/[ 18 0] = 
560 (Wilson & Rood 1994) are often used in a similar analysis (e.g., 
|Teyssier et al. 2002) . These values lead to the ratio [ 13 CO]/[C 18 0] 
7.3. 
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Fig. 4. Same as Fig.[3]but for the 1.3-km s" 1 component. The black contours in the upper panels, from left to right, are plotted as 
follows: from 0.7 (2.1 a) to 7.0 K km s" 1 (27 cr) in steps of 0.7 K km s" 1 ; at 0.5, 1, 1.5, and 2 km s" 1 ; and from 0.5 to 2.5 km s" 1 in 
steps of 0.5 km s" 1 . In the lower panels, these contours are from 0.3 (l.lcr) to 2.7 K km s" 1 (9.6<x) in steps of 0.3 K km s" 1 ; at 0.5, 
1, 1.5, and 2 km s" 1 ; and at 0.5, 1, and 1.5 km s" 1 . The C 18 peak position, from which the 1.3-km s" 1 component spectra were 
extracted (Sect. 3.4), is indicated by a yellow box in the left panels. The effective beam size 30" is shown in the lower left corner of 
the upper left panel. 




5 h 42 ffl 50' 40" 

RA 

Fig. 6. Same as Fig. [5] but across 7.7-10.41 km s" 1 of the 
C 18 0(2 - 1) line. The contours are as in Fig. [5] 



Once the optical thickness is determined, T QX can be calcu- 
lated using the radiative transfer equation [see, e.g., Eq. (1) in 




85,660 85.680 85.700 85.720 85.740 
Angle (cleg) 



Fig. 7. Position- velocity (PV) diagram of the 13 CO(2 - 1) emis- 
sion towards the environment of SMM 3. The PV slice runs from 
NW to SE, across SMM 3, as shown in Fig. [3] The black con- 
tours go from 10 to 90% of the maximum value of 3.6 K, in steps 
of 10%. The grey contours represent negative values. 



Paper I] . In this calculation, we assumed a beam filling factor of 
unity, and that the background temperature is equal to the cosmic 
background radiation temperature of 2.725 K. 

The derived values of T13, rig, and T QX are listed on Cols. (8) 
and (9) of Table [2] In this table, we use the symbol To for the 
optical thickness to denote its peak value. The quoted uncertain- 
ties were derived from the uncertainties of the peak intensities. 
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The 13 CO line appears to be optically thick, while the C 18 line 
shows a moderate optical thickness of rig - 0.8+q 3 . 

The total beam-averaged 13 CO and C 18 column densities 
were computed following Eq. (4) of Paper I. The spectroscopic 
parameters needed in the analysis, such as the electric dipole 
moments and rotational constants, were adopted from the JPL 
database. The derived column densities are listed in the last col- 
umn of Table [2l The quoted uncertainties were propagated from 
those associated with T ex , Av, and To (the average value of the 
i -errors of T QX and To were used). 




Velocity (km/s) 



C peak (east) 

13, 




Velocity (km/s) 

Fig. 8. The 13 CO(2 - 1) and C 18 0(2 - 1) spectra extracted to- 
wards SMM 3 (upper panel) and the position indicated in Fig.[4l 
(lower panel). Hyperfine- structure fits to the 13 CO lines, and sin- 
gle Gaussian fits to the C 18 lines are overlaid in green. The red 
dashed line overlaid on the spectra towards SMM 3 indicates 
the systemic velocity of SMM 3 as derived from C 17 0(2 - 1) in 
Paper III. In the lower panel, the red dashed line shows the radial 
velocity of the C 18 0(2 - 1) line. 



3.5. CO depletion in SMM 3 

In Paper III, we derived a CO depletion factor of fo = 10.8 + 2.2 
towards SMM 3 through C 17 0(2 - 1) observations (2778 re- 
solution). Another estimate of in SMM 3 can be obtained 
from the current C 18 data. We prefer to use C 18 for this anal- 
ysis rather than 13 CO, because C 18 emission is more optically 
thin than the 13 CO emission. Therefore, C 18 is expected to 
trace gas deeper into the core's envelope, and being less affected 
by foreground emission. 

The values of To and T ex for the C 17 0(2 - 1) line were de- 
rived to be 0.05 ± 0.03 and 11.0 ± 1.0 K, respectively (Paper 
III). As the oxygen -isotopic ratio [ 18 0]/[ 17 0] is only 3.52 
( |Frerking et al. 1982| ), it can be assumed that the C 18 0(2 - 1) 
line is also optically thin. Under the assumption of optically 
thin emission (r <?c 1), and adopting T QX - 11 ± 1 K, the 
C 18 column density computed from the integrated intensity is 
N(C 18 0) = 8.5 ± 1.7 x 10 14 cm" 2 . 



By smoothing the LAB OCA map to match the resolution 
of the C 18 map (30"), and regridding it onto the same pixel 
scale, the 870-yum peak flux density towards the (0", 0") po- 
sition of the C 18 map is determined to be 0.85 Jy beam -1 . 
Using the assumption that T& ust = T^ n , and that the dust 
opacity per unit dust mass a t 870 jim is /C87o = 1.7 cm 2 g" 1 
( |Ossenkopf & Henning 1 994; Paper I), we estimate that the cor- 
responding H2 column density is A/^tb) = 2. 5 + 0. 3x1 22 cm -2 
(see Papers I and III for further details). Therefore, the fractional 
abundance of C 18 towards the (0", 0") position is estimated to 
be jc(C 18 0) = N(C ls O)/N(U 2 ) = 3.4 ± 0.8 x 10" 8 . 

To compute fiy from C 18 data, we need an estimate of the 
"canonical", or undepleted, abundance of C 18 0. Using the stan- 
dard value 9.5 x 10" 5 for the abundance of the main CO isotopo- 
logue in the solar neighbourhood ( |Frerking et al. 1 982), we can 
write 

T 18 01 1 
4C 18 0) can = x(CO) can x^ = 9 .5xl0- 5 x— = 1.9X10" 7 . 

(3) 

The value of is then determined to be = 
x(C 18 0) can /x(C 18 0)obs = 5.6 ± 1.3. This is comparable 
within a factor of two to the value derived from C 17 0(2 - 1) 
data. The agreement is quite good given all the assumptions 
used in the analysis. We note that the (0", 0") position of the 
C 18 map is not exactly coincident with our previous C 17 
observation target position but the two are within the beam size 
of both the observations. 

Following the analysis presented in Miettinen (2012; 
Sect. 5.5 therein), the CO depletion timescale in SMM 3 is es- 
timated to be only Td ep ~ 1 - 3.7 x 10 4 yr [using the values 
r kin = 11.3 K and n(U 2 ) = 1.1 - 4.0 x 10 5 cm" 3 ; see Table[Q. 
This provides a lower limit to the age of the core. 



4. Discussion 

4.1. On the non-detection of molecular outflows 

Because SMM 3 is an early Class source, it is expected to 
drive a bipolar molecular outflow. The outflows can manifest 
themselves in broad non-Gaussian spectral-line wings. One of 
our original attempts of the present study was to search for out- 
flows driven by SMM 3 through 13 CO observations. However, 
in our data, there is no evidence for a large-scale 13 CO out- 
flow. In the Spitzer/IRAC 4.5-yum image of SMM 3 (Fig. [Q 
bottom panel), there are some 4.5-yum emission features that 
could be signatures of shock-excited material around SMM 3 
(e.g. JSmith & Rosen 2005llDe Buizer & Vacca 2010b . For com- 
parison, the Class 0/1 protosta r IRAS 05399-01 21 in Orion B9, 
which drives the HH 92 jet ( |Bally et al. 2002| ), exhibits spec- 
tacular IRAC 4.5-yum features along its parsec-scale bipolar jet. 
Clearly, higher resolution observations, and better outflow trac- 
ers, such as 12 CO and SiO, would be needed to clarify the out- 
flow activity of SMM 3. 



4.2. Low-velocity gas emission 

As was discussed in Paper II (Sect. 5.7 therein), the lower- 
velocity line emission seen towards Orion B9 is likely to come 
from the "low-velocity part" of Orion B, which probably ori- 
ginates from the feedback from the massive stars of the nearby 
Ori OB lb association. This fraction of the gas is likely to be 
located a few tens of parsecs closer to the Sun than the "regular" 
9-km s" 1 gas component ( Wilso n et al. 20051 ). 
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Table 2. Parameters of the spectra extracted from the selected positions. 



Position 


Offset 


Transition 


V LSR 


Av 




J T MB dv b 


to 




N 




") 




[km s- 1 ] 


[km s" 1 ] 


[K] 


[Kkms" 1 ] 




[K] 


[10 15 cm" 2 ] 


Main velocity component 




















SMM3 


(0,0) 


13 CO(2 - 1) 


8.11+0.09 


0.95 ± 0.36 


2.73 ± 0.45 


2.61+0.67 [7.17, 9.21] 












C 18 0(2 - 1) 


8.57 ± 0.05 


0.63 + 0.11 


2.12 + 0.25 


1.39 + 0.25 [8.00, 9.14] 


0.2±0.1 c 


11.0+ 1.0 C 


0.85 ± 0.2 d 


Lower-velocity component 




















C 18 peak 


(+85.5, +14.0) 


13 CO(2 - 1) 


1.52 + 0.06 


0.81+0.12 


5.17 + 0.61 


7.05 + 0.86 [0.21,2.38] 


6 3 +2Ae 

0*3 




18.2 + 6.8 






C 18 0(2- 1) 


1.38 + 0.06 


0.97 + 0.15 


2.75 + 0.61 


2.89 ± 0.47 [0.07, 2.52] 


9.5 ± 1.8^ 


2.8 + 1.0 



Notes. (a) Offset from the map centre in arcsec. ( ^ Integrated intensity is derived by integrating over the velocity range indicated in square 
brackets. (c) This optical thickness is estimated from that derived for C 17 0(2 - 1) in Paper III, and the reported T ex value is that of C 17 0(2 - 
1). ( ^ Derived from the integrated intensity under the assumption of optically thin emission. (e) Derived from the line intensity ratio, as described 
in Sect. 3AS f) Derived using the T MB and r values of the C 18 line. 



The C 18 column density towards the selected C 18 peak is 
2.8 ± 1.0 x 10 15 cm -2 . We can convert this to an estimate of the 
H2 column density as 



/V(H 2 ) = 



[H 2 



[ 12 CO] 1R 
x J x Af(C 18 0) . 



[ 12 CO] [C 18 0] 



(4) 



Using again the values [ 12 CO]/[C 18 0] = [ 16 0]/[ 18 0] = 500 and 
[ 12 CO]/[H 2 ] = 9.5 x 10~ 5 , we obtain N(U 2 ) = 1.5 ± 0.5 x 10 22 
cm -2 . This shows that the "low- velocity part" of Orion B also 
consists of dense gas, which conforms to our previous detection 
of, e.g., deuterated molecular species at comparable LSR veloc- 
ities. We also note that the prestellar core SMM 7 and Class 
protostar IRAS 05413-0104 seen towards Orion B9 have such 
low systemic velocities (~ 3.7 and -1.5 km s" 1 , respectively) 
that they are likely to be members of the low-velocity Orion 
B. Despite the estimated high column density of the clump-like 
13 CO/C 18 feature seen in the left panels of Fig. |4l it was not 
seen in LABOCA 870-yum emission at the ~ 3<x level (reflecting 
the strong CO depletion in the main velocity component). 

4.3. On the origin of the velocity gradient and implications for 
the core/star formation in Orion B9 

In Paper II, we speculated that the Orion B9 region has proba- 
bly been influenced by the feedback from the nearby Orion OB 
association, or more precisely, from the Ori OB lb subgroup. 
As mentioned above, this feedback is b elieved to be respon sible 
for the "low-velocity part" of Orion B (IWilson et al. 2005b . We 
believe that the discovery of a velocity gradient in the present 
study supports the possibility that Orion B9 region is affected by 
stellar feedback. 

To better illustrate the larger-scale view of the surroundings 
of Orion B9, in Fig. [9] we show a wide-field Herschel/SPIRE 
250-yum imag^Zl. A visual inspection of the image suggests 
that the Orion B9 region might be situated in quite a dy- 
namic environment. Towards the south, there is the active 
massive star-forming region NGC 2024 some 40' (~ 5.2 pc 
at 450 pc) from Orion B9. The northwest- southeast oriented 
250-yum filament emanating from NGC 2024 is likely related 
to the dense molecular ridge of NGC 2024 running in the 



7 Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA dPilbratt et al. 20101 . 
Orion B was observed as part of the "Hersch el Gould Belt Survey 
(GBS)" dAndre et al. 20101 using th e PACS ( |Poglitsch et al. 2010) 
and the SPIRE ([Griffin et al. 2010 ) instruments. For more de- 
tails, see http://gouldbelt-herschel.cea.fr The data 
are available from the Herschel Science Archive (HSA) at 
http : //herschel . esac . esa . int/Science_ Archive . shtml 



same direction (e.g., IThronson et al. 19841 IVisser et al~1 998; 
Wat anabe & Mitchell 20081) . The expanding H11 region of NGC 
2024 appears to be interacting with the molecular ridge 
dGaume et a l. 1992 ). The majority of the members of Ori OB lb 
association li e towards the west/ northwest from Orion B9 (see, 
e.g., Fig. 1 in lWilson et al. 2005b . 

It is apparent from the bottom panel of Fig. [9] that the Orion 
B9 cores, with the possible exception of SMM 7, belong to a 
common northeast- southwest oriented filamentary structure. As 
indicated in the figure, the cores at the northeastern part of the 
region have a lower LSR velocity or show multiple velocity 
components as found in our previous papers. As discussed ear- 
lier, a gradient of increasing radial velocity from NW to SE is 
seen across the mapped area, and SMM 3 appears to lie on the 
border of the velocity jump [cf. the case of the Class proto- 
star HH211 in Perseus/IC348 (ITobin et al. 201 1M . Interestingly, 
this velocity gradient appears to be orthogonal to the direc- 
tion of the 250-yum filament oriented NE-SW (SMM 3 also lies 
on the border of the dust filament). This raises the question 
whether this border could be tracing a shock layer of the inter- 
acting/colliding flows, within which there is a jump in the ve- 
locity of the gas. Such interaction might have been responsible 
for triggering the formation of SMM 3, and (some) of the other 
cores in the region. It seems more likely that, instead of the in- 
fluence of Ori OB lb, the feedback from the nearby expanding 
NGC 2 024 H11 region could have compressed the inital cloud 
region (Fukuda & Hanawa 2000|. Later, the cloud under a high 
pressure gradient may have fragmented into dense cores, out of 
which some, such as SMM 3, were collapsed into protostars se- 
quentially. 

SMM 3 showed the highest level of CO depletion among 
the cores studied in Paper III. This conforms to the fact 
that it also appears to be the densest core in Orion B9. For 
the core collapse induced by compression, the simulations by 
Hennebelle et al. (2003) suggest that the combined duration of 
the prestellar+Class phase is ~ 3.2 x 10 5 - 1.3 x 10 6 yr (de- 
pending on the rate of compression). As mentioned earlier, the 
estimated CO depletion timescale in SMM 3 is ~ 1 - 3.7 x 10 4 
yr, which is shorter than the core evolution timescales quoted 
above. On the other hand, the fragmentation timescale of the 
core is expected to be comparable to the signal crossing time, 
which for SMM 3 is estimated to be T cr0S s = D/cr 3D ~ 6.6 x 10 5 
yr (the projected core diameter across the 3.3cr contour is 86" 
or D = 0.19 pc, and 0-3D = V3<xid is the three-dimensional 
velocity dispersion). This agrees well with the above theoretical 
core lifetimes. From these considerations, we suggest that the 
formation of SMM 3, and of some other cores in Orion B9, was 
triggered by feedback from NGC 2024 (dynamical compression) 
some several times ~ 10 5 yr ago (cf. Fukuda & Hanawa 2000). 
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To better understand the velocity structure of the region on 
larger scales, larger maps of molecular-line emission would be 
needed. We note that the star formation in Orion B9, if trig- 
gered by stellar fee dback, mi ght resemble the situations in the 
p Qphiuchus (e.g., [Nutter et al. 2006b and the B59/Pipe Nebula 
JPeretto et al. 2012b , where the star formation appears to be in- 
duced by the feedback from the Scorpius OB association. 

5. Summary and conclusions 

A 4' x 4' region around the Class protostar SMM 3 in Orion 
B9 was mapped in 13 CO and C 18 7 = 2-1 lines with the 
APEX 12-m telescope. Our main results and conclusions can be 
summarised as follows: 

1. Both lines exhibit two well separated velocity components: 
one at ~ 1.3 km s" 1 and the other at ~ 8.7 km s" 1 . The latter 
is near the systemic velocity of SMM 3. The low- velocity 
component was already recognised in our previous studies, 
and it is believed to be related to the low-velocity part of 
Orion B. 

2. The 13 CO and C 18 emission are relatively widely dis- 
tributed compared to the dust continuum emission traced by 
LABOCA. The LABOCA 870-yum peak position of SMM 3 
is not coincident with any strong 13 CO or C 18 emission, 
which is in accordance with the high CO depletion factor de- 
rived earlier by us from C 17 0(2 - 1) (f D * 10.8). The CO 
depletion factor derived from C 18 data is within a factor 
of two from the previous estimate, i.e., fo - 5.6. No evi- 
dence for a large-scale outflow activity, i.e., high velocity 
line wings, was found towards SMM 3. 

3. The lower- velocity (~ 1.3 km s" 1 ) 13 CO and C 18 emission 
are concentrated into a clump-like feature at the eastern part 
of the map. We estimate that the H2 column density towards 
its C 18 maximum is ~ 10 22 cm -2 . Therefore, the lower- 
velocity gas seen along the line of sight is of high density, 
which is consistent with our earlier detection of, e.g., deute- 
rated molecular species (DCO + , N2D + ). 

4. We observe a velocity gradient across the 13 CO and C 18 
maps along the NW-SE direction (some hint of that is also 
visible in the lower- velocity line maps). Interestingly, SMM 
3 is projected almost exactly on the border of the velocity 
jump. The sharp velocity-gradient border provides a strong 
indication that it represents an interaction zone of flow mo- 
tions. 

5. We suggest a possible scenario in which the formation of 
SMM 3, and likely some of the other dense cores in Orion 
B9, was triggered by an expanding Hn region of NGC 
2024. This collect- and-collapse -type process might have 
been taken place some several times 10 5 yr ago. The NGC 
2024 region is known to be a potential site of induced , se- 
quential star formation (e.g., Fukuda & Hanawa 2000, and 
references therein). The case of Orion B9 suggests that we 
may be witnessing the most recent event of self-propagating 
star formation around NGC 2024. Larger- scale molecular- 
line maps would be needed for a better understanding of the 
larger- scale velocity structure of the region. 
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Fig. 9. Top: A wide-field Herschel/SPIRE 250-yum far-infrared image towards the Orion B9 and NGC 2024 (Flame Nebula) star- 
forming regions. The image is overlaid with black contours of LABOCA 870-yum dust continuum emission. The contours go from 
0.1 (~ 3.3<x) to 1.0 Jy beam -1 in steps of 0.1 Jy beam -1 . The green rectangle outlines the 4' x 4' area mapped in 13 CO(2 - 1) and 
C 18 0(2 - 1) in the present study. The white arrows indicate the cores which have a lower LSR velocity than the "main 9-km s" 1 
part" of Orion B9 or show multiple velocity components. The Class protostar IRAS 05413-0104 in the northeast corner was at the 
border of our LABOCA map and only weakly detected (see Fig. 1 in Paper I). Note how the Orion B9 cores are associated with a 
NE-SW oriented filamentary structure, and the cores with lower radial velocity (or with multiple velocity components) lie at the NE 
part of the structure. A scale bar indicating the 2 pc projected length is shown in the bottom left, with the assumption of a 450 pc 
line-of-sight distance. Bottom: A zoomed-in view of the upper image towards Orion B9. Selected cores are labelled. The LABOCA 
contours are as in the upper panel. 
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